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The interaction between a superconducting vortex in a type II superconducting film (SC) and a 
ferromagnet (FM) with out-of-plane magnetization is investigated theoretically within the London 
approximation. The dependence of the interaction energy on the FM-vortex distance, film thickness 
and different geometries of the magnetic structures: disk, annulus(ring), square and triangle are 
calculated. Analytic expressions and vectorplots of the current induced in the SC due to the presence 
of the FM are presented. For a FM disk with a cavity, we show that different local minima for the 
vortex position are possible, enabling the system to be suitable to act as a qubit. For FMs with 
sharp edges, like e.g. for squares and triangles, the vortex prefers to enter its equilibrium position 
along the corners of the magnet. 

PACS numbers: 74.78.-w, 74.25.Qt, 74.25.Ha. 

I. INTRODUCTION 

The interaction between superconductivity and magnetism has drawn a lot of attention in the last decades. To 
study the effects due to the interplay of the superconducting order parameter and the non-homogeneous magnetic 
field resulting from a ferromagnet (FM), several experimental groups fabricated periodic arrays of magnetic dots 
and antidots positioned above or under the superconducting filmi2ii4 Such ferromagnetic dots act as very effective 
trapping centers for the vortices which leads to an enhancement of the critical current. Recently, it was predicted 5 
that an increase of the pinning effects by two orders of magnitude can be realized in this way. After substantial 
progress in the preparation of regular magnetic arrays on superconductors and considering the importance of such 
structures for magnetic device and storage technologies, these hybrid systems became very interesting both from a 
theoretical and an experimental point of view. Macroscopic phenomena have already been explored experimentally, 
but a theoretical analysis of the magnetic and superconducting response in such systems is still in its infancy. 

In previously proposed models for the superconducting film (SC) interacting with a ferromagnet (FM) on top of 
ffe&LSiii the magnetic texture interacts with the SC current, which subsequently changes the magnetic field. The 
authors used the London approximation to describe this system since the sizes of all structures are much larger 
than the coherence length £. The thickness of the SC film and of the FM was assumed to be negligibly small (i.e. 
d « £,A). Elementary solutions for the interaction between the circular magnetic dot (bubble) or annulus (ring) 
with a vortex were found. Further, the creation of additional vortices by the FM near the surface was described®, 
by a simple comparison of the free energies of the system with and without the vortex. However, the spontaneous 
creation of a vortex-antivortex pair as a possible lower energy state was never considered. 

Other theoretical studies involving finite size ferromagnets were mainly restricted to the problem of a magnetic dot 
with out-of-plane magnetization embedded in a superconducting film^fl*ii Marmorkos et aZi£ were the first to solve 
the non-linear Ginzburg-Landau (GL) equation numerically, with appropriate boundary conditions for an infinitely 
long ferromagnetic cylinder penetrating the superconducting film, and found a correspondence between the value of 
the magnetization and the vorticity of the most energetically favorable giant-vortex state. The vortex structure of a 
SC disk with a smaller magnetic disk on top of it was numerically calculated in Ref^ using the first non-linear GL 
equation, i.e. neglecting the effect of the screening currents on the total magnetic field. Interesting vortex-antivortex 
k>( [ configurations and an interplay between the giant-vortex and multi-vortices were found. 

Most recently, the pinning of vortices by small magnetic particles was studied experimentalljji&Mii 5 . which was 
a motivation for our recent theoretical study of this system^. In the latter study we approximated the magnetic 
field profile by a magnetic dipole. In the present paper, we generalize these results in order to include the realistic 
magnetic field profile of the FM which in the present approach can be of arbitrary shape. The superconducting film 
lies in the xy-plane while the FM is positioned a distance I above the SC, and is magnetized in the positive z-direction 
(out-of-plane). To avoid the proximity effect and the exchange of electrons between the FM and the SC we assume a 
thin layer of insulating oxide between them as is usually the case in the experiment. 

The paper is organized as follows. In the next section we present the general formalism. In Sec. Ill, we discuss the 
pinning potential of the magnetic disk and magnetic annulus (ring) with out-of-plane magnetization. Further, the 
vortex-magnet interaction energy and supercurrent induced in the superconductor are calculated analytically and the 
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profiles are shown. We use these results in Sec. IV to investigate the manipulation of vortices in the case of a more 
complicated geometry of the magnet, i.e. magnetic disk with an off-center hole(s). In Sec. V the pinning properties of 
the magnet with square or triangular shape are analyzed and the most favorable trajectory of the vortex with respect 
to the magnet edge is determined. The influence of edges (corners) of the FM on the pinning is then discussed and 
our conclusions are given in Sec. VI. 



II. THEORETICAL FORMALISM 

We consider a ferromagnet of arbitrary shape with homogeneous out-of-plane magnetization M, placed outside a 
type II SC film interacting with a single vortex in the SC. Within the London approximation, the direct interaction 
energy between the vortex and the FM in a stationary magnet-superconductor system is given byl^ 

Umv = Y c j dVil) ' **] ~ \ j dVUm)f iv ■ M, (1) 

where $„ = (<5> p , $ z ) = (0,& /(2irp),0) denotes the vortex magnetic flux vector (<& is the flux quantum). The 
first integration is performed over the volume inside V^ 1 ' the superconductor, while V^ m%> in the second integral 
denotes the volume of the ferromagnet. Indices v and m refer to the vortex and the magnet, respectively, j denotes 
the current and h the magnetic field. 

The interaction energy in this system consists of two parts: (i) the interaction between the Meissner currents 
generated in the SC (j m ) by the FM and the vortex, and (ii) the interaction between the vortex magnetic field 
and the FM. In Ref^ we showed analytically that in the case of a point magnetic dipole (MD) for both in- and 
out-of-plane magnetization these two contributions are equal. Due to the superposition principle, the finite FMs 
with homogeneous magnetization can be represented as an infinite number of dipoles. Consequently, in our case of 
out-of-plane magnetized FM, the vortex-magnet interaction energy equals 

U mv = - f dV^h v -M. (2) 



In order to obtain the current induced in the superconductor by the ferromagnet, one should solve first the equation 
for the vector potential^ 

rot (rotA^ + -^9 (d/2 - \z\) A m = AnrotM. (3) 

This equation is rather complicated to be handled for a finite size FM, but, the analytic expressions for the induced 
SC current in an infinite superconducting film with thickness d (— | < z < |) in the MD case (with magnetic moment 
m) are knowni& 

■mdr ,_cm$ y m -y f f ( d\\ 2 



j™(x,y,z)= 2nX r R dqexpi-ql\z m \--)}q 2 J 1 (qR m )C(q,z), (4a) 



.,„„, <•///<!>,, •'' ■'•■„ /' , J / 1 d\> 2 



with 



j™ a (x,y,z)= 27rA3 " / dq e xpi-q\\z m \--)}q z J 1 (qR m )C(q,z), (4b) 



_ k cosh(fc(f + z)) + g sinh(fc(f + z)) 
[<1, Z) {k 2 + q 2 ) sinh(fcd) + 2kq cosh(fcd) ' 1 ' 



where k = yjl + q 2 , R m = y/ (x — x m ) 2 + (y — y m ) 2 is the distance between the dipole and the point of interest, and 
J v (a) is the Bessel function. The coordinates (x m ,y m , z m ) denote the position of the dipole. The magnetic moment 
of the magnet is measured in units of mo = &q\, and all distances are scaled in units of A. These units will be used 
in the rest of the paper. 

To find the supercurrent induced by a finite size FM above the superconductor, we make use of the superposition 
principle and consequently the above expressions 1 ll'ill have to be integrated over the volume of the ferromagnet. 
Thus, the value of the current is given by (a = x, y) 

j a (x,y,z)= [ C d (x,y,z)dV^ m \ (6) 
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III. MAGNETIC DISK (RING) - VORTEX INTERACTION 



In this section, we investigate the interaction between a vortex in an infinite type II superconducting film with 
thickness d (— ^ < z < i) and the magnetic disk with radius R and thickness D with out-of-plane magnetization, i.e. 

M — Q(R — p)0(zi — z)Q(z — zo)me^/(R 2 nD) (in units of M — <&/\ 2 ) located at distance / above (under) the SC 
(z =l,zi = l + D). 

Inserting the well known expression for the magnetic field of a vortex outside the SG-£ 

, / x £$0 Idqq , . / 2\z\-d\ 

^' Z ' = 2^A2 J ~Q~ ( q ' 6XP \ q 2 J ' ^ ' 

o 



L$ 



h vp (p, z) = -^—^ sgn(z) / -pp Ji (gi?) exp ( -q 



dq q 



Q 



2\z\-d 



(7b) 



into Eq. we find the expression for the magnetic disk-vortex interaction 

MR<5>1 



A 



Ux (Pv) , 



where p — p v denotes the position of the vortex, and 

oo 

U ± (p v ) = J dq^Jt (qR) J (q Pv ) E(q, I, D), 



(8a) 



(8b) 



where Q = k(k + qcoth(kd/2)), and E (q, I, D) = e" 9 ' ( e - qD - l) . 

For the case of a thin ferromagnetic disk above the thin superconducting film (d, D < 1), the following asymptotics 
can be obtained: 

(1) for p v < R, we found 



U±_ ( Pv ) 



dDR 



1 



1 



l + Vl 2 + R 2 4 (p+R2f/ 2 J 1 



(9a) 



(2) for p v > R and p v ~ I, 



U±_ ( Pv ) 



dDR 



R 2 2l 2 - p 2 v 



VP + 7v 8 (i 2 + P 2 v f 2 , 



(9b) 



(3) for p v > R,l, 



U ± ( Pv ) 



dDR 



1 



dn 

Vi 2 + pl~~ 



Pvd 
2 



p v d 
~~2~ 



(9c) 



Here, H v {x) and Y v (x) denote the Struve and Bessel function, respectively. 

In the case of a FM on top of the SC (I = 0), the above asymptotics l(Ha-c)) can be expressed in a more precise 
way: 

(1) for p v < R, 



dD^fp 2 v 



(2) for p v > R, 



Uj_ (p v 



dDR f 1 



7T {Pv 
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rv rp I 
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(10a) 



(10b) 
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K(x) and E{x) are the complete elliptic integrals of the first and second kind, respectively. Further expansion of the 
asymptotic behavior of the energy at large distances (Eqs. 10;) and dS>)) gi yes U± (p v ) ~ —DR/dp\. 

When we take the derivative of the interaction energy, Eq. ijHJl, over p v we obtain the force acting on a vortex in 
the presence of a magnetic disk: 



with 



^2 -L 



(11a) 



OO 

F±(p v ) = J dq^J 1 (qR)J 1 (qp v )E(q,l,D). 



(lib) 



For the case of a thin FM on top of a thin SC (d, D < 1, / = 0), we derived the following asymptotics 
(1) for p v < R, 



(2) for p v > R, 



F± (p v ) 



dD 
~ttR 



J? I \~ dD 
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irp v 
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Pvd 
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(12a) 



(12b) 



The latter expression reduces in the extreme p v >> R limit to F± (p v ) ~ —SDR/dp^, which is consistent with the 
asymptotic behavior of the interaction energy. 

The results for the full numerical calculation of Eq. (JSJ are shown in Fig. 1(a) for a magnetic disk with radius 
R = 3.0, and three values of the thickness D — 0.1, 0.5, 1.0, fixed total magnetic moment m — 1.0 (M = m/V^" 1 ^), 
and at distance I — 0.1 above the SC with thickness d = 0.1. The energy is expressed in units of Uq = $>q/tt\) and 
the force in F = ^o/ 7 ^ 2 - The 

magnet- vortex interaction energy increases if the magnet is made thinner, since the 
magnetization in that case increases. Also the magnetic field of the disk becomes more peaked near the magnet edge. 
In Fig. 1(b), the dependence of the interaction energy on the thickness of the SC is shown. The increased thickness of 
the SC makes the interaction stronger, due to the stronger response of the SC to the presence of the magnet. Notice 
that increasing the thickness beyond d >> A does not influence the energy (dashed curve in Fig. 1(b), Q w k(k + q) 
in Eq. ©). 

The vortex is attracted by the magnetic disk when the magnetization and the vortex are oriented parallel, inde- 
pendently of the value of the parameters. The interaction energy has its minimum just under the center of the disk, 
which is the energetically most favorable position of the vortex. The force acting on the vortex is purely attractive 
and it has its maximum at the edge of the magnetic disk (see Fig. 1(c)). For large distances between the magnetic disk 
and the vortex the interaction approaches the value obtained earlier for the case of a magnetic dipole»i& Note that in 
the limit R — > and D — > 0, Eq. JHJ corresponds to the out-of-plane dipole case of Ref>i&. In Refs&S the interaction 
energy between a magnetic nanostructure and a vortex in a thin superconductor [d « 1) was calculated. In their 
case, the thickness of the magnet was not taken into account (assumed to be infinitely thin), the superconducting film 
was taken very thin (d « A) and the FM was placed on top of the SC (in the same z — plane) which corresponds 
to I = in our case. In these limits, our equations reduce to the same expressions for the interaction energy like 
those given in Refsi&?. But note that the analytical expressions are still not completely reached in the distance 
range shown in Fig. 1 and therefore a numerical calculation of the full integral is necessary in order to obtain the 
magnet-vortex interaction energy. Therefore, from this point of view, our expression offers much more information 
(non-zero thickness of both FM and SC, and arbitrary position of the FM above the SC) without a real increase of 
the complexity of the numerical calculation. 

To better understand the attractive magnet-vortex behavior in this system, we calculated the supercurrent induced 
in the SC due to the presence of the magnet. As explained in Sec. II, this current can be obtained after integration 



FIG. 1: The magnetic disk-vortex interaction energy as function of the distance between the vortex and the center of the FM 
disk: for several values of (a) the thickness of the FM disk, and (b) the thickness of the SC. (c) Plot of the FM-vortex force, 
and (d) a vectorplot of the current induced in the SC due to the presence of the FM (same parameters as in (c)). The grey 
semicircle in (d) indicates the position of the edge of the FM. 
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of Eq. (JSJ . In the case of a flat magnetic disk, it has only an azimuthal component and reads 

oo 

UM = S ^T 1 J dqJ 1 (qR)J 1 (qp)E(q,l,D)C(q,z), (13) 


where p — y 1 x 1 + y 2 , and C(q,z) is given by Eq. J5J. For a FM placed under the SC, one should replace z by — z. 
The vectorplot of the current is shown in Fig. 1(d). One should notice that the direction of the current is the one 
normally associated with an antivortex (the clockwise direction) and that the current is maximal at the magnetic 
disk edge. This agrees with our previous results: the direction of the current explains the attraction between the FM 
disk and the vortex, and the position of the maximum of the current corresponds to the maximal attractive force. 
The problem is cylindrically symmetric, and a vortex approaching the magnet from any direction will be attracted 
for parallel alignment and repelled in the anti-parallel case. This important point was not fully explained in Ref£. 

Using the same procedure, for a magnetic annulus (Ri < p < R a ) with thickness D and out-of-plane magnetization 
(inset in Fig. 2(b)) we have 

M = ^—&(p - R t )®(Ro - p)Q{zx - z)Q(z - z Q )mr z , (14) 

' ann 

resulting in the vortex-magnet interaction energy 

U mv = J dq—J (qp v ) (RoJ^qRo) - RiJ^qRi)) E(q,l,D). (15) 

o 

The interaction energy and force curves for the magnetic annulus-vortex interaction are given in Fig. 2(a,b) and 
are in qualitative (but not quantitative) agreement with earlier results of Refi, which were obtained in the limit of 
an extremely thin SC and FM, namely d,D << A. Please notice again that in our calculation finite thicknesses of 
both SC and FM are fully taken into account. The most important result is that, in this case, the annulus-vortex 
interaction energy has a ring-like minimum, under the magnet. The exact radial position of this minimum depends on 
the SC parameters, the thickness of the magnet and its distance from the SC. The force acting on the vortex shows 
dual behavior - attractive outside the equilibrium ring and repulsive inside. 



FIG. 2: The out-of-plane magnetized annulus-vortex interaction: (a) the interaction energy, and (b) plot of the FM-vortex 
force. The contourplot of the interaction energy is shown as inset in (a) (dark color illustrates low energy, as will be the case 
throughout the paper, the dashed white semicircles illustrate the edges of the magnet). A schematic outline of the magnetic 
annulus is shown in the inset of (b). 

Due to the dual behavior of the FM-vortex force one naively expects different current flow in the superconductor 
inside and outside the annulus. We use again Eq. 10 and for the current induced in the superconductor we obtain 

oo 

3<p(p> z ) = ^T L / dq{R Jx{qR ) - RiJx{qRi)) Ji{qp)E{q,l,D)C {q,z) . (16) 
o 

In Fig. 3 we show a comparison between the currents induced in the SC in the case of a magnetic disk (dashed curves) 
and a magnetic annulus (solid curves), for different SC-FM vertical distances. When the magnet is positioned far 
above the superconductor, the vortex qualitatively does not feel the difference between the disk and the ring case, 
and the current induced in the SC shows a similar behavior (Fig. 3(a)). When approaching the superconductor, the 
influence of the central hole in the ring becomes more pronounced (Fig. 3(b,c)), and eventually the current changes 
sign (Fig. 3(d)). 

Obviously, the qualitative behaviors of all quantities outside the annulus approach those for the case of a magnetic 
disk. However, inside the ring, the situation is different. The nature of the magnet-vortex force changes and while 
the current flows in the clockwise direction outside the ring, inside the superconductor the direction of the current is 



FIG. 3: Comparison between the current induced in a SC (in units of jo = c<l?o/7rA 3 ) by a magnetic disk (dashed curves) and 
a magnetic annulus (solid curves) with the same outer radius, as the FM-SC vertical distance / decreases (a)-(d). 
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anti-clockwise in the case of a small FM-SC distance (i.e. l/X — 0.1, Fig. 3(d)). Please notice that due to the fact 
that the finite thickness of the SC is included in our calculations, the SC current contains also a z-dependence (i.e. 
through the C(q,z) function). 

From a look at Eqs. 1|15[) and Q16[) . one can see that the problem of a magnetic annulus actually can be modelled 
by two concentric magnetic disks with different radius and opposite magnetization. The problem is linear, and this 
will facilitate the calculation in the cases of non-cylindrically symmetric FMs. 

IV. MANIPULATION OF THE EQUILIBRIUM VORTEX POSITION WITH A MAGNETIC DISK 

CONTAINING A CAVITY 

In previous section, we discussed the pinning of vortices by a magnetic disk or annulus (ring) . We showed that the 
most energetically favorable position for the vortex is under the center of the magnetic disk (for parallel alignment) 
or under the annulus (equilibrium ring). Here we generalize the latter system and displace the hole in the disk from 
its central location. 

Referring to the previous section, we may consider this problem as a superposition of effects of two magnetic disks 
with opposite magnetization. The smaller radius magnetic disk with the opposite magnetization models the hole in 
the larger disk. The parameters of the magnet are: the outer radius R , the radius of the hole Ri, the center of the 
hole is at (xh,Uh) and the thickness of the FM is denoted by D. Therefore, using Eqs. I|8I13[) for two magnetic disks, 
one with radius R Q , and the other with radius Ri, with opposite magnetization and centered at (x,y) = (0,0) and 
(zhiUh), respectively, we investigate the pinning properties of such FM. 

The results of this calculation are shown in Fig. 4 for a magnetic disk with R a = 3.0, Ri — 1.0, (x^, yh) — (0.5, 0), 
/ = 0.1 and D = d = 0.5. In Fig. 4(a) we show the contourplot of the FM-vortex attractive force. It is clear that there 
are two local energy minima along the y = direction, where the force equals zero: in front of and behind the hole. 
In Fig. 4(b) the plot of the interaction energy is given along this direction for three positions of the hole Xh — 0.5, 0.8, 
and 1.0. The important result is that the two minima are not equal: the one closer to the outer edge of the magnet 
has higher energy (metastable state) and the one near the magnet center is the actual ground state. However, due to 

FIG. 4: The interaction of the vortex with a magnetic disk with an off-center hole (see left figure for schematics of the magnet 
configuration): (a) contourplot of the interaction energy for R — 3.0, Ri = 1.0, (xh,Vh) = (0.5,0), I = 0.1 and D = d = 0.5, 
with M — Mo (dashed circles indicate the edges of the magnet), and (b) the energy along the :r-axis, for three different locations 
of the off-center hole in the magnetic disk. 

the presence of the hole, the equilibrium position of a vortex is not exactly in the center, and depends on the position 
of the hole. The magnet is not cylindrically symmetric and we have two separate energy minima instead of a ring of 
minima as in the case of the magnetic annulus. Also, one could argue that a slowly moving vortex in a system with no 
temperature fluctuations could be trapped at the metastable position. Anyhow, the hole in a magnetic disk appears 
to be a powerful tool for a possible manipulation of the vortex position. However, one question arises: since there are 
two minima present in the interaction energy, is it possible to have two equilibrium states with the same energy? 

In order to construct such a situation, we introduced a second hole in the magnet, at a symmetrical position to the 
first one with respect to the center of the magnet. As an example, we took the parameters of the magnet as R = 3.0, 
Ri(for both holes) — 0.5, (xh,yh) = (±1.0,0). The interaction energy along the x-axis is given in Fig. 5. Two equal 
minima near the outer edge of the disk are found next to the magnet holes. However, the global minimum is still 
under the center of the magnet (see inset of Fig. 5) . 

To eliminate this minimum we allow the holes to touch each other and to form an "eight-hole" in the center of 
the magnet. The resulting interaction energy (Fig. 6(a)) has now only two equal minima along the y = direction, 
outside the hole, near the magnet edge, and a plateau-like behavior in the central region. Still, these minima are not 
the lowest energy states. The central global minimum of the interaction energy from the previous case is now split by 
the joined holes into two minima along the x = direction (see Fig. 6(b)). The latter minima represent the ground 
state for a vortex in the presence of a magnetic disk with an "eight-hole" . The vortex has two absolutely equal ground 
states and the same probability of eventually sitting in one of those. Thanks to this feature, a possible use of this 



FIG. 5: The magnetic disk with two symmetrical holes: plot of the FM-vortex interaction energy illustrating the position of 
the metastable and ground vortex states with respect to the position of the holes (dashed vertical lines) . The parameters of 
the system are R a = 3.0, Ri(for both holes) = 0.5, (x h ,y h ) = (±1.0,0), I = 0.1, D = d = 0.5 and M = M . 
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FIG. 6: The interaction energy of a vortex with a magnetic disk containing an "eight-hole": (a) interaction energy for the 
vortex positioned along the :r-axis at the y = line, (b) along the y-axis for x = 0, and (c) contourplot of the FM-vortex 
interaction energy (dark color- low energy) . In (c) the vectorplot of the current induced in the SC in the presence of the FM is 
superimposed. The solid circles denote the position of the holes in the FM above the SC. 



system for quantum computing can be analyzed, similarly to the quantum systems proposed before (see, for example, 
Ref.ll). 

In Fig. 6(c) the contourplot of the magnet- vortex interaction energy is given, together with a vectorplot of the current 
induced in the SC. The position of the eight-hole is denoted by the thick solid circles. Around the magnetic disk, the 
SC current flows in a clockwise direction, illustrating the general attraction between the FM and the vortex. However, 
under the magnet, current shows a dual behavior, and a vortex-antivortex-like current flow can be seen. Namely, at 
the equilibrium vortex-states we find "antivortex" current profiles while under the holes of the magnet a vortex-like 
current motion is present. This suggests the possibility that such a magnetic field configuration, for sufficiently strong 
magnetization, could induce interesting vortex-antivortex configurations if placed near a superconducting film. 

Using this approach, interaction of a vortex with magnets of more complicated shapes can be investigated. We have 
shown that the magnetic disk with a cavity is a nice example of how to control the vortex by the magnet-geometry. 
From the point of view of practical vortex manipulation, it would be interesting to move the vortex by changing the 
parameters of the system. Helseth^ proposed a system in which a magnetic disk (with magnetization Mi) is placed in 
the center of a magnetic ring (with magnetization M 2 ), where the disk and the ring can have opposite magnetization. 
In Fig. 7 we show the outline of the system (upper inset), and the calculated interaction energy and force acting on 
the vortex, for the parameters W\ — 0.5A e , W2 — 0.75A e and W3 = 1.0A e and D = 0.1A on top of a SC with thickness 
d = 0.1A. Here A e = A 2 /d denotes the effective penetration depth. We suppose a thin oxide layer between the SC and 
the FM with thickness 0.01A. First we consider the Mi = M2 case. Helseth claimed that in this case, a slowly moving 
vortex will be attracted and sit under the annulus. The plot of the interaction energy in Fig. 7(a) shows two energy 
minima, namely, under the disk, and under the annulus (ring-energy minimum, see inset of Fig. 7(a)). From our 

FIG. 7: Disk-ring magnetic structure (see inset of (a) for a schematic view of the configuration) above the SC: (a) FM-vortex 
interaction energy for M\ = M2, (b) for M2 = —Mi. The contourplots of the energy are given as insets (dashed lines indicate 
the edges of the magnet). The force acting on the vortex is shown in (c), as the solid curve for case (a) and dashed one for case 
(b). The parameters are Wi = 0.5, W2 = 0.75, W3 = 1.0, D = 1.0, and I = 0.1, all in units of the effective penetration depth 
A e = X 2 /d. 

calculations it is clear that the ground state for the vortex would be in the center of the magnetic structure and not 
under the annulus as claimed in Ref£. However, one could argue that a vortex, slowly moving towards the magnet, 
could rest in a metastable state under the ring, if there are no fluctuations in the system. In Fig. 7(c), the solid line 
shows the force acting on the vortex. Please note that in our calculation, the positions of the extremes correspond to 
the magnet edges, while in Rcf. 9 this was not the case. Also, the peaks in the forces in our calculation are finite, due 
to the finite thickness of the magnet. It should be noted that the relation between the two energy minima and the 
acting forces strongly depends on the parameters and it can change in favor of the annulus if it is made wider. 

In Rcf. 9 it is also stated that in the opposite case, when the annulus magnetization changes sign, i.e. M2 = —Mi, the 
vortex is attracted to the center of the disk. The interaction energy we calculated in this case is plotted in Fig. 7(b) 
(the force is given in Fig. 7(c) by the dashed line). It is clear that the annulus forms an energetic barrier which 
prevents a "slowly moving" vortex from reaching the central position. However, the global minimum of the energy is 
under the center of this magnetic structure (the energy shown in Fig. 7(b) decreases monotonously and reaches zero 
at infinity), implying that the vortex would definitely sit under the magnet for the considered parameters. Different 
values for R\, R2, R3 could make the central minimum higher, and the energy barrier would then be able to repel the 
vortex. To conclude, in order to use this magnetic structure for vortex manipulation, one should not only overcome 
the experimental difficulties to realize such a structure, but also be careful about the influence of the parameters on 
the behavior of this system. 

V. SQUARE AND TRIANGULAR MAGNETIC DISKS: INFLUENCE OF THE CORNERS ON THE 

VORTEX PINNING 

Up to this point, we have only considered the interaction of a vortex with magnets having circular symmetry, 
namely magnetic disks, rings, and combinations of those. Here we consider FMs with broken circular symmetry, 
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namely square or triangular magnetic disks. We put the center of our Cartesian coordinates in the center of the 
magnet and the x axis along one of its sides. In this case, Eq. can not be solved analytically, and triple numerical 
integrations must be performed. In the case of a rectangular FM we have 

2 oo A/2 B/2 

Umv{x v ,Vv) = J ^ E(q,l,D) J dx J dy Jo (qRo) , (17) 

-A/2 -B/2 



where R v = y/Jx — x v ) 2 + (y — y v ) 2 - A and B are the dimensions of the magnet (with thickness D) in the x and y 
direction, respectively. As before, the distance between the FM and the SC surface is I. 

The components of the current induced in the SC are obtained in an analogous way, by numerical integration of 
Eq. |J?jJ, using the expressions I|4al5|l . 

In Fig. 8(a) the interaction energy with a vortex is shown, along two directions (lines of symmetry): (i) diagonal 
(dotted line); (ii) horizontal (solid line). In the inset, the contourplot of the energy is given. For comparison, we 
also give the results for a magnetic disk (thick dashed curve) with R 2 tt = AB and the same thickness D. As far 
as the pinning potential of a square magnetic disk is concerned, the asymmetry is rather small and the result is not 
much different from the one of an equivalent circular disk. Only in the region near the edge of the magnet, some 
discrepancy between the pinning potential of the square magnet in the diagonal direction and the disk approximation 
is observed. Moving further from the magnet, this discrepancy disappears. The energetically favorable position of 
the vortex is under the center of the magnet. The situation is somewhat different for rectangular-shaped FM disks. 
In Fig. 8(b) we show the results for A = AB = 2.0 and, in this case, the corresponding circular disk becomes a very 
poor approximation. Far away from the magnet, this approximation becomes better, as expected. 



FIG. 8: Rectangular magnetic disk above a SC film: (a) FM-vortex interaction energy for a square magnet case, namely with 
sides A = B = 1.0A, thickness D — 0.1A and magnetization M = Mo, at / = 0.1 above the SC with thickness d = 0.1A (the inset 
shows the contourplot of the energy, white line indicates the edge of the magnet), and (b) for A = 4B = 2.0. For comparison, 
the FM-vortex interaction energy in the case of the magnetic disk with the same volume is given by the thick dashed curves. 

From the behavior of the interaction energy we have seen that the vortex is attracted to the center of the square 
or triangular magnet for parallel orientation of the magnetization and the vortex magnetic field. This corresponds 
qualitatively to the case of a magnetic disk. However, the broken circular symmetry of the magnet introduces some 
changes in the magnet-vortex interaction. In Fig. 9 we show the contourplot of the force acting on the vortex and the 
vectorplot of the current induced in the SC, both for the case of a square and triangular magnetic disks. It is obvious 
that the attractive force is stronger at the sides than at the corners of the magnet. Therefore, the vortex approaching 
the magnet at the side of the magnet will be attracted stronger than on the diagonal direction. As far as the current 
is concerned, it has the direction associated with an antivortex. Near the magnet, the current follows the shape of the 
magnet and is maximal along the sides of the magnet. Further from the square or triangular magnet, the behavior 
of the current is more similar to the case of the circular magnetic disk. In Fig. 9(c) the rectangular magnet case is 
shown. One important feature should be noticed: the FM-vortex attraction force is stronger at the longer side of the 
rectangle. 

As one can see in Fig. 9, the maxima of the FM-vortex interaction force are located on the sides of the magnet. 
Therefore, one may expect that the energetically preferable direction of vortex motion when attracted by the FM 
is perpendicular to its sides. To investigate this, we put the vortex in different positions (open dots in Fig. 10(a)) 
and follow its trajectory using molecular dynamics simulations. The results are shown in Fig. 10(a). If the initial 
position of the vortex is along the lines of symmetry of the magnet, the vortex moves straight towards the center of 
the magnet. Otherwise, the trajectory of the vortex is distorted towards the corner of the magnet. It actually appears 
that the vortex avoids the maxima of the attractive force. This is counterintuitive but can be explained by the profile 
of the FM-vortex interaction energy, given in Fig. 10(b,c). One should notice the "wave" shape of the energy (if going 
along the ring around the magnet, see fig. 10(b)). Following a circle around the magnet, the energy has its minima at 
the corners of the magnet (denoted by black triangles in Fig. 10(b,c)) and the saddle points are on the sides (white 
triangles). The periodicity in Fig. 10(c) corresponds to the number of corners of the ferromagnetic disk. 



FIG. 9: The contourplots of the FM-vortex attractive force for parallel orientation of the magnetization and the vortex magnetic 
field, for: (a) an equilateral triangular magnetic disk, with sides a = 1.0A, (b) square magnet, with sides A = B = 1.0A. The 
other parameters correspond to the ones in Fig. 8. The vectorplots of the current are superimposed. In (c), the contourplot of 
the force for the rectangular magnet case is given, A — 1.6A, B — 0.6A. Dark color represents high force intensity. Positions of 
the edges of the magnets are given by white solid lines. 
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From Fig. 10 it is obvious that the interaction energy for any position of the vortex lowers steep towards the 
center of the magnet, but also towards the corners. This induces the distortion of the vortex trajectory and gives the 
impression that the vortex approaching the magnet from the corners is more energetically favorable. 

FIG. 10: The trajectory of the vortex when attracted by the triangular magnetic disk (same parameters as in Fig. 9(a)): 
(a) vortex paths with respect to the attractive force landscape (the edge of the magnet is illustrated by white line), (b) the 
contourplot of the triangular magnet-vortex interaction energy (dark color - low energy), and (c) the interaction energy along 
the ring indicated by dashed line in (b). 



VI. CONCLUSION 

We applied the London theory to investigate flux pinning in SC films due to the presence of a ferromagnet situated 
above (or under) the SC, where the finite thickness of both FM and the SC were fully taken into account. In the case 
of a magnetic disk or annulus (ring) , we obtained analytic expressions for the FM- vortex interaction energy, force and 
the screening currents. We also derived the asymptotic behavior of the interaction potential and the force for specific 
values of the involved parameters. In the case of a magnetic disk with an off-center hole we showed the existence of 
two local minima in the FM-vortex interaction energy - the ground state and the metastable one. By changing the 
position of the hole, the position of the equilibrium moves with respect to the magnetic disk center. We also showed 
that in the case of a FM disk with two touching holes ("eight-hole") two minima with equal energy but different 
vortex position appeared. The probability of a vortex sitting in one of these two states is the same, which makes this 
system interesting as a possible qubit. To further investigate the influence of the magnet geometry on its pinning 
properties, we calculated the pinning potential for square and triangular shaped ferromagnets. A substantial breaking 
of the circular symmetry occurs and the attractive force acting on the vortex is stronger at the sides of the magnet 
than at the corners. Also, making one side of the rectangular magnet longer enlarges the attractive force along it with 
respect to the other side. Although counterintuitive, we showed that the vortex approaches the non-circular magnet 
rather along the corners than perpendicular to the sides, following the gradient of the potential. 
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